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ABSTRACT 

Observations with the Wide-Field Camera of the Hubble Space Telescope (HST) are presented for eight 
intrinsically luminous quasars with redshifts between 0.16 and 0.29. These observations, when combined with 
a similar HST study of the quasar PKS 2349-014, show that luminous nearby quasars exist in a variety of 
environments. 

Seven companion galaxies brighter than M v — — 16.5 ( H 0 = 100 km s Mpc , D 0 = 1.0) lie within a pro- 
jected distance of 25 kpc of the quasars; three of the companions are located closer than 3" (6 kpc projected 
distance) from the quasars, well within the volume that would be enclosed by a typical L* host galaxy. The 
observed association of quasars and companion galaxies is statistically significant and may be an important 
element in the luminous-quasar phenomenon. 

Apparent host galaxies are detected for three of the quasars: PG 1116 + 215, 3C 273, and PG 1444 + 407; 
the hosts have an average absolute magnitude of about 0.6 mag brighter than L* . The agreement between the 
previously published major-axis directions in ground-based images and in the present HST images of 3C 273 
and PG 1444 + 407 constitutes important evidence supporting the reality of these candidate host galaxies. 

Upper limits are placed on the visual-band brightnesses of representative galactic hosts for all the quasars. 
These limits are established by placing galaxy images obtained with HST underneath the quasars and measur- 
ing at what faintness level the known galaxies are detected. On average, the HST spirals would have been 
detected if they were as faint as 1 mag below L*, and the early-type galaxies could have been detected down to 
a brightness level of about L*, where L* is the Schechter characteristic luminosity of field galaxies. Smooth, 
featureless galaxy models (exponential disks or de Vaucouleurs profiles) are fitted to the residual light after a 
best-fitting point source is subtracted from the quasar images. The results show that smooth spiral galaxies 
brighter than, on average, about L*, would have been detected. These upper limits, or possible detections, are 
consistent with, for example, the eight luminous quasars studied in this paper, occurring in host galaxies that 
have a Shechter luminosity function with a lower cutoff of in the range 0.01-0. 1L* . 

Tests are performed to determine if our failure to detect, in some cases, luminous host galaxies could be an 
artifact caused by our analysis procedures. These tests include comparing the measured point-spread function 
(PSF) for our HST observations with the PSFs used in previous ground-based studies of host galaxies, mea- 
suring the fluctuations in the sky signals that were subtracted from the quasar images, evaluating empirically 
the effects of using different stellar PSFs in the analysis, carrying out the subtraction of the stellar (nuclear) 
source in different ways, creating and analyzing artificial active galactic nuclei (AGNs) with known surface 
brightnesses, and fitting the observed quasar light to an analytic model that includes a host galaxy. Our 
analysis procedures successfully pass all these tests. 

Subject headings: galaxies: clusters: general — quasars: general 


I, INTRODUCTION 

It has long been been recognized that quasars might reside 
in galaxies and that the quasar phenomenon might be strongly 
influenced by the presence of companion galaxies. Depending 
upon your theoretical bias, you may expect to find quasars in 
old galaxies, in young galaxies, or in interacting galaxies, all 
possibilities that have been advocated in the published liter- 
ature. Very close companions might provide fuel for, or initiate 
processes that lead to, the quasar phenomenon. 

1 Based on observations with the NASA/ESA Hubble Space Telescope , 
obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. 

2 Also Institute for Advanced Study, School of Natural Sciences, Princeton, 
NJ 08540. 


Over the past two decades there have been a large number of 
ground-based observational programs in which the properties 
of host galaxies of quasars were investigated. In fact, a con- 
sensus view has been developed that quasars reside in lumi- 
nous galaxies, with the intrinsically bright quasars residing in 
the most luminous galaxies. Moreover, many authors have 
concluded that radio-bright quasars are in elliptical galaxies 
and radio-faint quasars are in spiral galaxies. Some representa- 
tive papers providing evidence for these views are Kristian 
(1973), Wyckoff et al. (1980), Wyckoff, Wehinger, & Gehren 
(1981), Tyson, Baum, & Kreidl (1982), Hutchings et al. (1982), 
Gehren et al. (1984), Heckman et al. (1984), Malkan (1984), 
Malkan, Margon, & Chanan (1984). Boroson, Persson, & Oke 
(1985), Smith et al. (1986), Hutchings (1987), Stockton & 
MacKenty (1987), Yee (1987), Hutchings, Janson, & Neff 
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(1989), Romanishin & Hintzen (1989), Veron-Cetty and 
Woltjer (1990), Hutchings & Neff (1992), Dunlop et al. (1993), 
and McLeod & Rieke (1994a, b). The analyses presented in 
these pioneering studies are difficult because of atmospheric 
seeing; the light from the bright central (nuclear) sources may 
be a few magnitudes brighter than the total emission from the 
host galaxies. 

In this paper, we concentrate on analyzing carefully the 
HST data and on providing new quantitative results. We gen- 
erally avoid speculating on whether our data are consistent 
with previous ground-based observations. We believe that the 
observers who made the ground-based observations are best 
able to decide on the extent of consistency or lack of consis- 
tency between our results and theirs and to make new obser- 
vations to test our results and theirs using the information 
provided by the HST images. We will do everything we can to 
facilitate the comparison of our results with previous and 
future ground-based results. 

Some of the papers that have previously discussed the com- 
panions of low-redshift quasars include Bahcall, Schmidt, & 
Gunn (1969), Bahcall & Bahcall (1970), Gunn (1971), Robinson 
& Wampler (1972), Burbidge & O’Dell (1973), Stockton (1978), 
Stockton (1982), Hutchings et al. (1982), Heckman et al. (1984), 
Malkan et al. (1984), Green & Yee (1984), and Yee (1987). The 
most dramatic results presented in the present paper refer to 
the very close galaxy companions that are not visible in pre- 
viously obtained ground-based images. 

The excellent optical characteristics of the repaired HST 
make possible improved observational studies of the hosts of 
small-redshift quasars. The characterization of faint galactic 
material (the host galaxy) in the vicinity of a bright stellar 
source (the quasar) can be done with increased reliability using 
the HST to minimize the possibility that quasar light is mis- 
taken for diffuse galactic light. Close galactic companions 
(within 2" or 3" of the quasar nucleus) may also be most con- 
vincingly detected with HST. 

We describe in this paper the systematic results that have 
been obtained for the first eight objects in a sample of 20 of the 
intrinsically most luminous ( M v < —22.9, for H 0 = 100 km 
s~ 1 Mpc~\ Q 0 = 10) nearby (z < 0.30) radio-quiet and radio- 
loud quasars selected from the Veron-Cetty & Veron (1993) 
catalog. Images of the first four of the quasars (PG 0953 + 414, 
PG 1116 + 215, PG 1202 + 281 [GQ Com], and PG 
1307+085) were presented in Bahcall, Kirhakos, & Schneider 
(1994, hereafter Paper I); images of four additional quasars 
(3C 273, PKS 1302- 102, PG 1444 + 407, 3C 323.1) have not 


previously been discussed (hereafter, “ the new quasars ’ ). 
Hutchings et al. (1994a) reported on observations of two active 
galactic nuclei (AGNs) imaged with the planetary camera (PC) 
of WFPC2. In the absence of an empirical stellar point-spread 
function (PSF), they used a model PSF and image restoration 
to analyze their images and to compare with ground-based 
data. 

Our principal results for all eight quasars and the hosts are 
summarized in Table 1, which lists the following quantities for 
each object: the date observed, the redshift, the distance in 
kiloparsecs that corresponds to an angular separation of 1" as 
seen from Earth, the apparent K-magnitude (from Veron-Cetty 
& Veron 1993), the absolute F-magnitude, and the range of 
galaxy magnitudes at which eight representative galaxies 
(selected from HST WFC2 observations with the same filter) 
could be detected when artificially placed underneath the 
quasar images. The last column gives, for all but three cases, 
the average limiting absolute magnitude at which the eight 
representative galaxies could have been detected as hosts. For 
the three quasars PG 1116 + 215, 3C 273, and PG 1444 + 407, 
we think it is likely that we have detected host galaxies; the last 
column of Table 1 contains for these three objects the host 
galaxy magnitudes obtained by fitting a theoretical model of 
an exponential disk to the residual light (see discussion in § 5). 
Accurate (1") coordinates and multicolor CCD photometry for 
six of the eight quasars are given in Kirhakos et al. (1994). 
Three of the quasars, 3C 273, PKS 1302— 102, and 3C 323.1, 
are radio loud; the other five quasars are radio quiet 
(Kellermann et al. 1989). 

The eight quasars studied in this paper have an average 
absolute visual magnitude of <M V > = -23.9, approximately 2 
mag brighter than the brightest galaxies of rich clusters 
(Hoessel & Schneider 1985; Postman & Lauer 1995) and about 
3.4 mag brighter than the characteristic (Schechter-) magnitude 
for field galaxies (Schechter 1976; Kirshner et al. 1983; Efsta- 
thiou, Ellis, & Peterson 1988). The average apparent brightness 
for the quasars is <F> = 15.2. At the typical redshift of the 
quasars, (z ~ 0.2), an If galaxy would have an apparent magni- 
tude of V « 18.6. 

The observations reported in this paper are all taken with 
the F606W filter of the WFC2 of the HST (see § 2 for a dis- 
cussion of the observations and relevant references). Since we 
do not know the shape of the spectra for the host galaxies of 
luminous quasars, we have chosen to neglect the relatively 
small difference between the V band and the F606W filter and 
also to omit all K -corrections. If the host galaxies were ordi- 


TABLE 1 


Sample of Completed Cycle 4 Quasars and Detectable Hosts 


Object 

Date 

(1994) 

z 

Distance“ b 

(kpc) 

V 

M K (QSO) b 

m » m ( host ) 

F606W 

<M, im (host)> b 

F606W 

PG 0953 + 414 

Feb 3 

0.239 

2.4 

15.3 

-24.1 

£ 18.4-20.0 

-20.0 

PG 1116 + 215 

Feb 8 

0.177 

1.9 

15.0 

-23.7 

£17.9-20.0 

— 20. 7 C 

PG 1202 + 281 

Feb 8 

0.165 

1.8 

15.6 

-23.0 

>18.6-20.3 

-18.9 

3C 273 

Jun 5 

0.158 

1.8 

12.8 

-25.7 

£16.4-18.8 

— 21.7 C 

PKS 1302-102 

Jun 9 

0.286 

2.7 

15.2 

-24.6 

£17.9-21.2 

-19.8 

PG 1307 + 085 

Apr 4 

0.155 

1.8 

15.3 

-23.1 

£17.7-20.5 

-19.3 

PG 1444 + 407 

Jun 27 

0.267 

2.6 

15.7 

-24.0 

£18.7-21.3 

— 20.8 C 

3C 323.1 

Jun 10 

0.266 

2.6 

16.7 

-22.9 

£18.6-21.3 

-18.6 


a Values correspond to an angular separation of 1" as seen from Earth. 

b Computed for Q 0 = 1.0 and H 0 = 100 km s 1 Mpc l . In this cosmology, brightest cluster galaxies have M v - 
- 22.0 and the characteristic (Schechter-) magnitude for field galaxies is MJ! = - 20.5. 


c Candidate host galaxy (see § 5). 
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nary elliptical galaxies, then the K-band magnitudes would be 
approximately 0.26 mag fainter at z = 0.20 (0.21 mag or 0.33 
mag fainter at z = 0.15 or z = 0.25, respectively) than the HST 
F606W magnitudes. For an Scd galaxy, the offsets would be 
somewhat smaller, 0.10 mag fainter at z = 0.20 (0.08 mag or 
0.13 mag fainter at z=0.15 or z = 0.25, respectively). We 
assume here that V = m( F606W). Thus, if the host galaxies of 
luminous quasars have spectra like normal galaxies (see 
Coleman, Wu, & Weedman 1980; B. McLeod 1995), then the 
actual K-magnitudes (and the actual K-absolute magnitudes) 
should be between 0.1 mag and 0.3 mag fainter than we quote. 

Are the HST observations described here sensitive enough 
to have detected luminous host galaxies? The answer is yes. 
We already have — in the same observational program — 
detected and measured with ease a faint nebulosity in the PKS 
2349-014 system which extends over about 50 kpc and has an 
absolute magnitude of only about M v = — 20.7 (see Bahcall, 
Kirhakos, & Schneider 1995). To put the question in per- 
spective, we note that a typical brightest cluster elliptical 
(M v = —22.0) at a redshift of z = 0.2 would produce about 
1 x 10 6 counts (detected electrons) in an annulus between 0'.'5 
and 5'.'0, i.e., between 1 and 10 kpc from the quasar center of 
light. The background noise from the sky and from read noise 
would only be about 10 3 counts in the same annulus, corre- 
sponding to a huge signal-to-noise ratio. The primary limi- 
tations on our detection of low surface brightness emission is 
not from the Poisson noise in the sky or from the read noise, 
but rather from imperfect sky subtractions, flat-fielding, and 
scattered light (extended PSF). These limitations are discussed 
in § 8; they are not severe enough to prevent our detecting L* 
and brighter host galaxies (for example, the luminosity 
observed in the PKS 2349 — 014 system; Bahcall et al. 1995). 

The quasar 3C 273 is a special case. It has a brighter appar- 
ent magnitude (V = 12.8, a brighter absolute magnitude 
(M y = -25.7), and a larger (~ 70 kpc total extent) and bright- 
er candidate host [Af K (host) ~ —21.7] than any of the other 
quasars in our current sample. It also has a unique optical 
counterpart for its well-developed radio jet. The host galaxy 
for 3C 273 is best seen in Figure 7. 

This paper is organized as follows. In § 2, we describe the 
observational procedures, including the determination of the 
stellar PSF. We outline, in § 3, the method by which we set 
limits on the brightnesses of representative HST- observed host 
galaxies. This section describes how we subtract the stellar 
PSF (§ 3.1), displays a catalog of observed HST galaxies that 
are used as simulated host galaxies (§ 3.2), and discusses the 
determination of the faintest magnitudes at which we could 
detect host galaxies underneath each observed quasar (§ 3.3). 
We discuss in § 4 fits of the residual diffuse light to smooth, 
featureless galaxy models, either an exponential disk or a de 
Vaucouleurs profile. We present in § 5 the characteristics of the 
candidate host galaxies that we have found for PG 1116 + 215, 
3C 273, and PG 1444 + 407. The HST images reveal separate 
galaxies whose projected separations from the centers of light 
of the quasar images are very small; these companion galaxies 
are discussed in § 6, and the properties are summarized in 
Table 6. McLeod & Rieke (1994b) have reported H-band 
detections of host galaxies for all eight of the quasars we 
discuss in this paper and have given sufficient detail that we 
can simulate the host galaxies they describe. We compare, in 
§ 7, our limits on the K-band brightnesses of simulated 
McLeod-Rieke galaxies with the reported H- band galaxy mag- 
nitudes and infer a minimum V - H color for the host galaxies 


if both observational analyses, theirs and ours, are correct. In 
§ 8, we attempt to answer the question: What have we done 
wrong? Why have we failed to detect the luminous host gal- 
axies reported by previous investigators who used data 
obtained with ground-based telescopes? We discuss in § 8 the 
difference between the PSF measured for our HST images and 
the PSFs used in previous ground-based studies of host gal- 
axies (§ 8.1), the sky subtraction with the HST images (§ 8.2), 
the results of experiments in which the images of different stars 
(observed at different positions in the camera or through differ- 
ent filters) were subtracted from each other (§ 8.3), a compari- 
son of separate analyses performed when the subtraction of a 
stellar point source was optimized on the diffraction spikes 
alone or on an azimuthally average light intensity (§ 8.4), and 
an analysis of simulated AGNs (§ 8.5). We summarize and 
discuss our principal conclusions in § 9. 

2. OBSERVATIONS 

Most of the details of the observation procedures are given 
in Paper I; we summarize here the main features of the obser- 
vations. We describe in § 2.1 the main features of the quasar 
images and present in § 2.2 the empirically determined stellar 
PSF. 

2. 1 . Quasar Images with F606 W 

The quasars were observed with the Wide-Field Camera 2 
(WFC2) through the F606W filter (see Burrows 1994a for a 
description of the instrument and the response curve of the 
filter). The F606W filter is similar to the V bandpass but is 
slightly redder; the mean wavelength and FWHM of the 
F606W system response are 5940 A and 1500 A, respectively. 
This filter was chosen because of its high throughput, which 
was the primary consideration for the present exploratory 
program. At a given angular radius, the scattered light in the 
WFC2 is about 5 times less than the scattered light in the 
planetary camera WFPC2 (see Krist & Burrows 1994). We 
chose to use the WFC instead of the PC in the original formu- 
lation of this program because of the likelihood that the host 
galaxies might have low surface brightnesses that extended 
over areas large compared to the WFC resolution (0?1 or 
about 0.2 kpc). 

The Space Telescope Science Institute (STScI) photometric 
calibration was verified by comparing, in quasar fields 
obtained at Palomar by Kirhakos et al. (1994), HST galaxy 
magnitudes with ground-based photometry of the same gal- 
axies. For the purposes of this paper, the V and the F606W 
photometric bands are sufficiently similar (Bahcall et al. 1994; 
Holtzman et al. 1995) that we will sometimes use V and 
m F606W interchangeably. 

All quasar images were located at a distance of 5" + 1'.'5 
from the center of Wide-Field CCD 3 and, except for PG 
1307 + 085, were observed for three separate exposures of 1 100 
s, 600 s, and 100 s; the exposures for PG 1307 + 085 lasted 1400 
s, 500 s, and 200 s. The image scale of the Wide-Field detectors 
is 0'.'0966 pixel -1 (Trauger et al. 1994), which means that the 
data are severely undersampled. The innermost regions of the 
quasar images are saturated in all of the exposures out to a 
radius %0''3. Typically, there were about nine saturated pixels 
in the shortest exposures and 50 saturated pixels in the longest 
exposures. The exposures for 3C 273 contained about 200 satu- 
rated pixels in the longest exposure and about 30 saturated 
pixels in the shortest exposure. 
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The initial data processing (bias frame removal and flat-field 
calibration) was performed at the STScI with their standard 
software package. The individual images of each quasar were 
aligned to better than 0.3 pixels; this made it easy to identify 
and eliminate cosmic-ray events. The flat-field corrections were 
based upon preflight calibrations; these calibrations remove 
the small-scale (few pixels) sensitivity variations. However, 
there are apparent large-scale sensitivity variations present in 
the calibrated data. The typical signal and rms of the noise of 
the sky in the long exposures (in detected photons per pixel) 
are 125 and 13, respectively. The sky level corresponds to a 
surface brightness of approximately 22.5 mag arcsec -2 . The 
observed noise is slightly larger than expected from a com- 
bination of shot noise and CCD readout noise (seven electrons 
per pixel). The formal detection limits for extended objects, 
calculated from these numbers, are extremely faint. However, 
one limiting factor for the detection of large-scale galaxies is 
the imperfect match between the data and the preflight flat- 
field calibrations. 

Figure 1 (Plates 13-14) displays the images of the eight 
quasars, with the only additional processing beyond that per- 
formed by STScI being the removal of cosmic rays. No lumi- 
nous host galaxies are visible in the raw data, which contrasts 
with what was expected based upon our studies of simulated 
AGNs with bright host galaxies (see § 8.5). 

2.2. The Stellar Point-Spread Function 

The point-spread function (PSF) was determined using a set 
of images of an isolated V = 10.5 mag star (F 141) in M67 that 
was placed on the corners of a box 3" on a side. The box was 
centered at a point displaced 2'.' 5 away from the central pixel of 
Wide-Field CCD 3 in each of the two orthogonal directions. At 
each corner of the box, a 12 s exposure was taken. In addition, 
a 2.6 s and a 70 s exposure were obtained at one corner of the 
box. The M67 observations, which were obtained on 1994 June 
5 (the same day the observations of 3C 273 were taken), cover 
the saturation range present in the quasar images. In these 
calibration images, the shortest exposures were saturated in 
the inner 0"2, the intermediate exposures were typically satu- 
rated to about 0'.'3, and the longest exposures were saturated 
out to 0"6. 

We did not obtain PSF calibrations at the same time as each 
quasar observation, since — given our limited observing time 
this would have reduced our sample size by a factor of 2. 
Simulations made prior to launch of what galaxy hosts HST 
could detect (based upon previously reported ground-based 
observations) indicated that the expected host galaxies would 
be easily detected and PSF subtraction would not be crucial. 
(The validity of these simulations was verified at least in part 
by the HST observations of PKS 2349-014; see Bahcall et al. 
1995). Fortunately, we have found (see the discussion in § 8.3) 
that the variation with time of the PSF does not significantly 
influence our inferred sensitivity limits. 

Figures 2 (Plate 15) and 3 show the point-spread function of 
the HST produced by the F606W filter; this PSF was con- 
structed by combining the short, intermediate, and long ex- 
posures of M67-F141. Figure 2 displays the two-dimensional 
empirical PSF generated from the three exposures of different 
lengths of the star F141 in M67. The two-dimensional profile is 
slightly asymmetric; the radial spokes are probably due to 
scattering within the instrument (see Krist & Burrows 1994). 
This empirical PSF has a larger dynamic range and contains 
much more signal in the wings than the images of the quasars. 



Fig. 3. — The azimuthally averaged F606W PSF. The figure compares the 
azimuthally averaged Burrows (calculated) point spread function ( solid line ) 
and the azimuthal average of the empirical point spread function shown in the 
previous figure ( dashed line). 

Figure 3 shows the azimuthal average of the empirical PSF 
for F141 and the F606W PSF calculated by Burrows (1994b). 
We measured the change in the profile caused by excluding the 
spikes; it is small, about 15%. 

The extended low surface brightness that is apparent in 
Figures 2 and 3 constitute the principal difficulty in investigat- 
ing faint host galaxies of quasars. Even with the diffraction- 
limited PSF of HST, if one fails to subtract precisely the stellar 
PSF from a quasar image, the residual stellar light from the 
imperfect subtraction may simulate the appearance of a host 
galaxy. The actual PSF at a fixed position in the camera is 
somewhat time dependent, since the focus changes due to 
desorption and to breathing modes of the telescope system. We 
will discuss in § 8.3 the limits on the detectability of faint hosts 
that result from imperfect PSF subtraction. 

The PSF determined using F141 is not optimal for removing 
the nuclear light from quasars, since FI 41 is considerably 
redder than the low-redshift quasars. Gilliland et al. (1991) 
measured the following colors for F141: V = 10.49, 

B—V = +1.11, U-B = +1.02, and V-R = +0.58 (see also 
Eggan & Sandage 1964; Chevalier & Ilovaisky 1991). The 
Veron-Cetty & Veron (1993) catalog gives accurate colors for 
five of the eight quasars discussed in this paper. The average 
quasar colors are B— V = +0.05 and LJ — B — —0.86. Since 
the measured PSF depends upon color, the empirical PSF 
constructed using images of F141 is not optimal for subtrac- 
ting nuclear light from quasars. Empirically, however, the color 
dependence of the stellar PSF does not cause large errors in the 
subtraction process (see § 8.3). 

3. BRIGHTNESS LIMITS ON ORDINARY HOST GALAXIES 

For the five cases in which we do not detect candidate host 
galaxies, we established the brightness upper limits that are 
given in Table 1 on host galaxies by artificially placing a series 
of observed and simulated galaxies underneath the quasar 
images. We then determined, with both objectives and subjec- 
tive (visual inspection of the data) techniques, the surface 
brightness levels at which the galaxies are detectable. The 
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detailed procedure used to establish a limiting brightness for 
each quasar-galaxy pair is described in this section. 

The first stage of the analysis, described in § 3.1, consists of 
subtracting a stellar PSF from the observed quasar image. We 
then construct a catalog, presented in § 3.2, of observed HST 
images that span the morphological types of host galaxies. The 
concluding step, described in § 3.3, is the termination of the 
limiting magnitude at which each of the putative host galaxies 
in our catalog of HST images could be detected underneath 
each of the eight observed quasars. 

3. 1 . Subtracting the Stellar PSF 

The best fit of the stellar PSF to the quasar image was 
determined with the downhill simplex method (described in 
Press et al. 1986) using the x 2 calculated in relatively narrow 
regions centered on the diffraction spikes. There are in prin- 
ciple three input parameters to the fit: the central location of 
the quasar (two coordinates) and the total brightness of the 
PSF. However, the location of the quasar could be determined 
accurately by eye (typically to 0.5 pixels, i.e., 0705) by adjusting 
the positions of the residuals from the diffraction spikes. In 
practice, we performed subtractions in which the central posi- 
tion was located by eye and, independently, was located by the 
simplex method; the two procedures gave indistinguishable 
results for the residual diffuse emission. Pixels that were satu- 
rated in either the quasar data or the PSF data, or that showed 
a large disagreement between the PSF and the quasar data (3 a 
from the mean deviation), were eliminated from the x 2 fitting. 

Since the template PSF contained a much higher signal in 
the wings than any of the quasar data, the PSF was always 
scaled down in brightness in subtracting the PSF from a 
quasar image. For all practical purposes, the PSF subtraction 
was effectively noiseless. 

Figure 4 (Plates 16-17) shows the result of subtracting the 
composite best-fit stellar image, Figure 2 (of M67 calibration 
star F141) from the long-exposure images of each of the eight 
quasars. Each panel in Figure 4 is 20" x 20". The contrast has 
been set to emphasize low surface brightness features; in the 
long exposures, the quasar images are saturated to a radius of 
»076. 

For the images shown in Figure 4, the intensity of the sub- 
tracted stellar images was adjusted visually so as to best elimi- 
nate the diffraction spikes. In an independent procedure, 
azimuthal averages of the residual intensity differences between 
the quasar and the stellar PSF were minimized using a x 2 
routine. The two methods gave quantitatively similar results 
(see § 8.4). We also explore in § 8.3 the accuracy of the PSF 
subtraction by subtracting different stellar images from each 
other. 

3.2. Images of T est Galaxies 

In order to quantify our detection limits for faint nebulosity, 
we established a catalog of representative images of individual 
galaxies that were added to the HST quasar images. Except for 
special tests, we used images of galaxies that were observed 
with the Wide-Field CCDs of the WFPC2 through the F606W 
filter. We note that there is a possible bias in using only the 
kinds of galaxies that were detected by HST (none of which are 
perfectly smooth) to determine our detection limits. For this 
reason, we also used smooth, simulated test galaxies (both 
exponential disks and de Vaucouleurs spheroids). Unless 
explicitly specified otherwise, all the results in this paper refer 
to the catalog of HST- observed galaxy images. 


Figure 5 (Plate 18) shows the collection of eight observed 
galaxy images that were placed underneath the quasars and 
whose surface brightnesses were scaled in order to determine 
empirically our sensitivity to extended nebulosity centered on 
the quasar. The first six images appear in Figure 2 of Paper I; 
they are reproduced here for the reader’s convenience. 

All eight of the galaxies in Figure 5 were taken from HST 
observations in this program; they represent a variety of 
galaxy types. Figure 5a represents a pair of interacting spiral 
galaxies, Figure 5b is an edge-on spiral, Figure 5c is a ringed 
spiral seen face on, Figure 5 d is a smooth elliptical, Figure 5e is 
an SO or an Sa, and Figure 5/ appears to be a disturbed barred 
spiral. The last two panels, Figure 5 g and Figure 5/i, contain 
pairs of interacting galaxies. 

The exponential scale lengths of the three spiral galaxies, 
Figures 5b, 5c, and 5/ are, respectively, 076, 076, and 078. We 
also used in the tests described below images expanded by a 
linear factor of three for the two spiral galaxies shown in 
Figures 5b and 5 f Therefore, the range of exponential scale 
lengths for the galaxies considered here are, at the redshifts of 
the quasars, between l and 5 kpc. The exponential scale length 
of the Galaxy is about 3.5 kpc (Bahcall 1986), a value which 
also describes well the average exponential scale length of 
other local galaxies and disk -dominated galaxies out to red- 
shifts of order 0.8 (see Mutz et al. 1994). 

The de Vaucouleurs effective radii for the three elliptical 
galaxies shown in Figures 5d-5e are, respectively, 372, 474, and 
172. At the redshifts of the quasars considered in this paper, the 
effective radii vary between 2 and 9 kpc. For the Galaxy, the 
effective radius of the spheroid is about 2.7 kpc (Bahcall 1986). 
It is difficult to define uniquely a length scale for the interacting 
galaxies shown in Figures 5g-5h; the major axes of these 
systems are respectively, 9 " and 6", corresponding to distances 
of order 12-18 kpc. 

Table 2 summarizes the characteristic length scales for the 
eight test galaxies shown in Figure 5. The exponential scale 
length are given in arcseconds and, for z = 0.2, in kiloparsecs. 

We have measured from other HST images the angular sizes 
of galaxies in a cluster at redshift 0.17 and verified that the 
cluster galaxies have similar angular sizes to the galaxies that 
were chosen to appear in Figure 5 and which were used in our 
sensitivity tests (see below). 

3.3. The Faintest Detectable Test Galaxies 

To determine the limiting apparent magnitude at which a 
specific test galaxy could be detected, the specified galaxy was 

TABLE 2 

Characteristic Length Scales 
of the Test Galaxies 


Galaxy 

Scale Length* 

5a 

0’.'6 b (1.3) b 

5b 

ore" (i.2) b 

5c 

3'.'2 C (6.7) c 

5d 

4"4 C (9.2) c 

5e 

1 72 c (2.5) c 

5f 

0'.'8 b (1.7) b 

5g 

r.'7 b (3.6) b 

5h 

273 b (4.9) b 


a Numbers in parentheses are 
scale lengths in kiloparsecs for 
z - 0.2. 

b Exponential scale lengths. 
c de Vaucouleurs effective radius. 
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Fig. 6. —The typical dependence of x 2 vs. apparent magnitude of a putative 
host galaxy. In the particular case illustrated here, the SO/Sa galaxy shown in 
Fig. 5e was added to the quasar PG 0953 + 414. 


added to each quasar image before subtraction of a stellar 
PSF. The galaxy sizes were held fixed as they appear in Figure 
5. In all the cases considered in this paper, the center of the 
putative host galaxy and the quasar were aligned. Flost gal- 
axies are more easily detected if they are not centered on the 
quasar (see Paper I). 

The x 2 for the best PSF fit was calculated as a function of 
galaxy brightness. As the galaxy image was made progressively 
fainter, the x 2 dropped until the galaxy signal became lower 
than the background noise. The template galaxies were taken 
from images that also contained the quasars and so have the 
same noise and PSF as the quasar images. The noise in the 
template galaxies were not scaled separately, since the original 
observations were in all cases rather close to the detection limit 
(in the worst within 50% of the original brightness). When the 
computed x 2 equaled twice the asymptotic value, a clear break 
point in the x 2 versus brightness curves, we could easily detect 
the galaxy placed underneath the quasar. Visual inspection of 
the subtracted images showed that in some cases (in particular 
for the edge-on spiral, Fig. 56), a galaxy could be detected % 1 
mag fainter than indicated by the x 2 test described above. 

Figure 6 shows the typical dependence that is found for x 2 as 
a function of apparent magnitude of the putative host galaxy. 


Table 3 gives, for all eight of the quasars analyzed, the limit- 
ing magnitude down to which each of the eight galaxies shown 
in Figure 5 could be detected when the galaxy was centered 
on the quasar image. The limiting magnitudes presented in 
Table 3 take into account the information obtained by visual 
inspection. 

For a given quasar, the limits of detectability of the galaxy 
images spanned a range of 1.5 + 0.3 mag (see Table 3). Table 1 
presents (in the next-to-last column) the magnitude ranges at 
which the faintest host galaxies, fainter versions of Figures 
5a-5Ji, could be detected. We averaged, for each quasar, the 
limiting magnitude to compute a representative limiting mag- 
nitude of a detectable host galaxy. The last column of Table 1 
gives for all eight quasars the average absolute magnitude at 
which the eight galaxies in Figures 5a-5h were visible when 
placed underneath that quasar. 

Averaging all 64 determinations of limiting detectability for 
the eight galaxies. Figures 5a-5h, placed underneath the eight 
quasars, we find that the ensemble-average apparent magni- 
tude at which host galaxies could be detected is 19.3, about 4.1 
mag fainter than the average quasar magnitude. The average 
absolute magnitude limit is 

^F606w( lim ) = - 19 6 . (0 

which is about 1 mag fainter than L* for field galaxies (see 
Table 1, note a). 

In general, spirals are easier to detect then ellipticals, 
because spirals have more features. On average, the spirals 
could be detected to a limiting magnitude of 19.4 mag, i.e., 

MF 606 wO im > spirals) = - 19.5 . (2) 

The two early-type galaxies, Figures 5d-5e, could only be seen 
down to 18.3, 

MF 606 wO in L early-type) = -20.6 , (3) 

1 mag brighter than the ensemble average. Typically, we can 
detect the interacting galaxies shown in Figures 5g-5h down to 
20.1, 0.8 mag fainter than the average of the other six galaxies 
shown in Figure 5. 

Some previous authors have suggested that the host galaxies 
of quasars were large in size. In order to simulate large host 
galaxies, we also studied scaled versions, 3 times larger but 
with the same surface brightness profiles as in Figure 5, of the 
galaxies in Figure 5b and Figure 5/ The results for the two 
simulated larger galaxies are given in the last two rows of 
Table 3. On average, we could detect these large putative host 
galaxies to 19.8 mag, about 0.5 mag fainter than the ensemble 
average of the eight observed galaxies shown in Figure 5. 


Galaxy 


5a .. 
5b .. 
5c... 
5d .. 
5e . . . 
5f ... 
5g .. 
5h .. 
5bx3 
5fx3 


0953 + 414 


19.7 
20.5 

18.7 

18.4 

18.8 

20.5 
20.1 
19.7 
19.4 
19.7 


1116 + 215 


19.1 

19.8 

18.8 

18.9 

17.9 

19.0 

20.0 
19.0 
19.3 
19.0 


TABLE 3 

Limiting Apparent Magnitudes 


1202 + 281 

3C 273 

1302-102 

1307 + 085 

1444 + 407 

3C 323.1 

19.7 

17.7 

19.4 

19.5 

20.3 

20.4 

20.1 

18.0 

20.0 

20.0 

20.4 

20.5 

18.8 

17.0 

18.7 

18.5 

19.7 

19.8 

18.8 

17.0 

18.4 

18.3 

19.6 

19.0 

18.6 

16.4 

17.9 

17.7 

18.7 

18.6 

20.0 

18.0 

19.5 

19.4 

20.4 

20.2 

20.3 

18.8 

21.2 

20.5 

21.3 

21.3 

19.8 

18.0 

20.2 

20.0 

20.0 

20.6 

19.5 

18.8 

20.2 

20.0 

20.5 

21.3 

19.7 

18.5 

19.9 

19.7 

20.8 

20.8 
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We are able to detect galaxies that are well-separated from 
the quasars down to about 24th magnitude, approximately 5 
mag fainter than the limiting magnitudes listed for host gal- 
axies in Table 3. 

4. FITS WITH MODEL GALAXIES 

Many authors have determined the magnitudes of host gal- 
axies by fitting the quasar images to a model consisting of a 
point source plus an analytic model of a featureless galaxy (an 
exponential disk or a de Vaucouleurs profile) (see, e.g., McLeod 
& Rieke 1994a, b; Dunlop et al. 1993; Veron-Cetty & Woltjer 
1990; and references therein). This procedure assumes that the 
host galaxies are smooth and featureless (which is different 
from most of the //ST-observed galaxies shown in Fig. 5 of 
§ 3), that the standard stellar PSF is a good match to the PSF 
of the nuclear region of the quasar, and that there are no bright 
companion galaxies. 

We shall see in this section that, not surprisingly, we obtain 
limits on, or possible measurements of, host galaxy magnitudes 
that are about 1 mag less stringent when we fit the residual 
light to smooth galaxy models than when we determine (as in 
§ 3) the limits of detectability of //ST-observed spiral galaxies. 
This difference also reflects the approximately 1 mag difference 
in the detectability limits found in the previous section (see eqs. 
[2] and [3]) between spirals and early-type galaxies. 

We have fitted combined models of a stellar nucleus plus a 
host galaxy with an exponential disk to each of the long- 
exposure quasar images. The quasar nucleus was approx- 
imated by the composite stellar PSF (Fig. 2). The general form 
of the assumed fit was 

Light(x, y) = C, x PSF(x, y) + C 2 

x Host(x, y; length scale) + Sky Offset , (4) 

where the PSF and the host galaxy intensity, Host, were nor- 
malized to unity when integrated over all positions, (x, y). We 
used analytic functions, an exponential disk or a de Vaucou- 
leurs profile, for the host galaxy surface brightness. 

For many, but not all, of the tests we ran, the central inten- 
sity and position of the stellar PSF were fixed so as to best 
eliminate the diffraction spikes (see Fig. 4). We varied the initial 
guesses for the amount of light in the stellar nucleus, C u over a 
wide range and also varied by large factors the initial guesses 
for the amount of light, C 2 , in the host galaxy. In addition, we 
varied the characteristic length scale, the effective radius (for a 
de Vaucouleurs profile) or the scale length (for an exponential 
disk), between 1" and 3". Moreover, we did experiments in 
which the fit was optimized within different annuli; typically, 
the annulus used was 075-478, but we also used other annuli 
(such as 179-7:7). 

The formal x 2 fits to the combined model of a smooth expo- 
nential disk plus stellar nucleus were generally poor; the 
typical x 2 per degree of freedom was 2-4, but in a number of 
cases the fit was much worse. Some of the mismatch was 
caused by the imperfect subtraction of the stellar PSF and is 
evident in Figure 4; we could not completely remove the dif- 
fraction spikes in the subtracted images of the quasars. In 
retrospect, it is perhaps not suprising that the residual light is 
not well described by a featureless, smooth analytic galaxy 
model. 

The results for the best-fitting apparent and absolute magni- 
tudes of the model host galaxies are shown in Table 4 for the 
special case in which the initial guess for the parameters was 


TABLE 4 

Total Magnitudes for a Model of a Point Source plus an 
Exponential Disk 


m F606W 

Quasar z (Exponential Disk) (Exponential Disk) 


PG 0953 + 414 0.239 19.5 -19.9 

PG 1116 + 215 0.177 18.1 -20.7 

PG 1202 + 281 0.165 18.7 -19.8 

3C 273 0.158 16.8 -21.7 

PKS 1302-102 0.286 18.4 -21.4 

PG 1307 + 085 0.155 19.0 -19.4 

PG 1444 + 407 0.267 18.8 -20.8 

3C 323.1 0.266 19.1 -20.5 


C i — C 2 — 0.0. In general, the fits with an exponential disk 
were rather stable to parameter variations, with a total range 
of only ±0.2 in the average absolute magnitude for the sample. 
The average absolute magnitude for the best-fitting exponen- 
tial disks of the host galaxies is 

F 606 wO* m )> 8 qiiasars = -20.5 ± 0.2 mag , (5) 

which is equal to the absolute magnitude of an L* galaxy. A 
similar range, ±0.2 mag, was found for the variations of most 
of the fits to the magnitudes of the assumed host galaxies in 
individual quasars, although occasionally differences as large 
as ±0.5 mag were found for individual model host galaxies. 

The magnitudes given in Table 4 and in equation (5) must be 
regarded with caution. We have no proof that the residual light 
being fitted by the smooth galaxy models is in the assumed 
form of a smooth exponential disk. In fact, the fits to this form 
are, as mentioned above, poor. Moreover, if we start the / 2 
optimization routine at a large negative value for the host 
galaxy intensity, i.e., C 2 is negative and large, we find host 
galaxies with negative fluxes. If we change the sign of the nega- 
tive fluxes obtained in this way, they correspond to host galaxy 
magnitudes of about the same average brightness as found 
when we made an initial positive guess for the light in the host 
galaxy, i.e., they are consistent with equation (5), The scale 
lengths for the cases with negative fluxes, however, are all very 
large, so large that the surface brightness of the galaxies is 
comparable to the noise in the sky. One indication of the size of 
the error in the fitting procedure may be the absolute magni- 
tudes obtained by inverting the sign of the negative fluxes, 
which is of order of M v = —20.5 mag. 

We also find it somewhat suspicious that all the exponential 
scale lengths cluster around the values 172 ± 073. This clus- 
tering suggests that the fits could be representing some instru- 
mental features instead of an astronomical property. 

The fits to a de Vaucouleurs profile were even less stable to 
variations of the initial parameters than were the fits to an 
exponential disk. The best-fitting effective radii were some- 
times very large, comparable to the outermost radius con- 
sidered in the fitting process. In a few of the cases in which a de 
Vaucouleurs profile was assumed, we found negative fluxes for 
the host galaxies even if the initial value of C 2 was taken to be 
zero. These cases had very small effective radii, typically of 
order 075, so that a large fraction of the light in the model host 
galaxy occurred in areas that were not part of the fitting region. 

We conclude that there is some evidence that there is 
residual light not associated with the stellar nuclei of the 
quasars that has an average absolute magnitude of about l*. 
However, this light is not well described by a smooth galaxy 
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model. In addition, changes in the initial guesses for the fitting 
parameters can sometimes produce negative galaxy fluxes. We 
conclude that our observations are not able to determine deci- 
sively whether smooth, featureless galaxies are present with 
visual luminosities about equal to L*. 

Observations with HST at different telescope roll angles are 
required in order to decide if the residual light that is fit in this 
section by smooth galaxy models is associated with the quasar 
or is somehow produced in the WFPC2. 

5. CHARACTERISTICS OF CANDIDATE HOST GALAXIES 

For the three quasars PG 1116 + 215, 3C 273, and PG 
1444 + 407, we have detected residual emission that is plausibly 
attributed to host galaxies. Tests described in this paper (see 
§ 8) and in Paper I show that we have detected diffuse emission 
in excess of the presumed stellar nuclear component (as 
modeled by the composite PSF shown in Fig. 2) for these three 
(and to a lesser extent for some of the other) quasars we have 
studied. However, with the exception of 3C 273 (for which 
there is a good agreement with previous ground-based obser- 
vations in the orientation of the host galaxy image), the avail- 
able HST observations are insufficient to establish definitively 
the origin of the light that is in excess of the nuclear emission. 
The light left over after the subtraction of a best-fitting stellar 
PSF could result from: (1) the faint host galaxies of quasars; (2) 
the difference between the composite PSF of the red star F141 
(Fig. 2) and the PSF of the blue quasar nucleus; or (3) excess 
light from companion galaxies that are projected very close to 
the quasars (see, e.g., the close companions of PG 1202 + 281, 
PKS 1302-102, and 3C 323.1 in Fig. 1 of Paper I and Fig. 4 of 
this paper). 

In this section, we present estimates of the characteristics of 
candidate host galaxies for PG 1116 + 215, 3C 273, and PG 
1444 + 407, based upon the assumption that all the detected 
residual light is contained in galactic hosts. We discuss 3C 273 
first, since the evidence is most convincing for this case. 

Figure 4 d, which displays for the bright radio quasar 3C 273 
the residual image after subtraction of a best-fit point source, 
shows evidence of an elliptical host. The galactic host is, 
however, most clearly seen when the image is compressed by 
combining the light in 3 x 3 pixels into a single pixel (see § 8.6 
for a more general discussion of this technique). 

Figure 7 (Plate 19) shows the clearer view of the host galaxy 
of 3C 273 that is apparent in the compressed image (in this 
case, without the subtraction of a stellar PSF). The candidate 
host galaxy is relatively bright, with an aperture magnitude 
(measured on the PSF-subtracted image) of m F606W « 17.1 
(M v = —21.4) in the annular region between 3" and 15" (i.e., 
between 5 and 27 kpc). The total measured extent of the host 
galaxy along the major axis is about 30" and about 20" along 
the minor axis. The agreement in the orientation of the candi- 
date HST elliptical host with what was inferred from previous 
ground-based images (see especially Tyson et al. 1982; and 
Roser & Meisenheimer 1991) of 3C 273 convinces us that this 
residual emission corresponds to a real host galaxy (see § 9). 

The isophotes of the residual nebulosity do not appear to be 
centered on the 3C 273 nuclear emission. The center of the 
quasar light is displaced about 2" (4 kpc) from the center of 
the brightest isophotes of the nebular emission and about 2" in 
the opposite direction from the center of the faintest iosphotes. 

It is difficult to determine the total luminosity of the host 
galaxy of 3C 273, since we are not able to measure the lumi- 
nosity in the innermost region. If (as has been done in analyz- 


ing the ground-based observations) we adopt a smooth, axially 
symmetric de Vaucouleurs profile centered on the quasar and 
fit to a one-dimensional average of the residual light (excluding 
the diffraction spikes and regions where the emission is 
saturated), we find V ~ 16.0 (M v = -22.5). However, the 
effective radius of the best-fit model is about 2", which implies 
that more than half the light in the model lies in the region in 
which no measurements can be made. 

For PG 1116 + 215 (see Fig. 46), it is difficult to decide what 
diffuse emission is due to the host and what is due to the 
extended PSF. There appears to be a faint (23.4 mag arcsec' 2 ), 
smooth ringlike structure with a radius of about 2" surround- 
ing almost half of PG 1116 + 215, plus an additional diffuse 
protrusion. If this nebulosity is real, it corresponds to a host 
galaxy with an aperture magnitude of m F606 w ~ }9.0 (0.1 mag 
brighter than the average limiting magnitude estimate for the 
eight HST-observed galaxies listed in Table 3), and a total 
diameter of ~ 20 kpc. 

We have performed aperture photometry for a number of 
different annuli surrounding the nucleus of PG 1116 + 215; the 
results are relatively robust. For an annulus with an inner 
radius of 2" and an outer radius of 7" (i.e., between 4 and 14 
kpc), we estimate a candidate host galaxy magnitude of 
m F 606 w « 18.6 (M v = -20.1). It is difficult to estimate how 
much of the light within 2" is associated with the candidate 
host galaxy. If one assumes that the surface brightness profile is 
flat within 2", then the additional light in this inner region is 
less than 0.2 mag. 

The candidate host galaxy associated with PG 1444 + 407 
can be seen clearly in the residual image of Figure 4. Inter- 
estingly, this nebulosity is not as obvious in the compressed 
images. The faint nebulosity that appears to surround PG 
1444 + 407 has an approximate surface brightness of 23.6 mag 
arcsec” 2 . The approximate magnitude of the nebulous feature, 
which may extend as far as 15 kpc, is m F606W % 19.4 (M v = 
— 20.3). The nebulosity is centered on the quasar. 

Table 5 lists the estimated characteristics for the three candi- 
date host galaxies. The columns contain the following informa- 
tion: the quasar name, the estimated aperture magnitude, the 
estimated magnitude determined by assuming that the host 
galaxy has an exponential disk profile, the absolute magnitude 
of the host galaxy computed from the aperture magnitude (or 
the exponential disk model) assuming that it has the same 
redshift as the quasar, and comments regarding the galaxy. 

We also measured the residual surface brightness in regions 
well separated from the quasar; the magnitudes we obtained 
by multiplying the average surface brightness by the approx- 
imate area over which the surface brightness extended agree 
well with the results of aperture photometry. This procedure 
yields a reasonable estimate for the diffuse emission only if the 
host galaxy profile is flat and if the residual light is all due to a 
host galaxy. The exponential disk magnitudes were determined 
by finding the parameters that best fit the quasar image 
assuming all the light is produced by the composite PSF of 
Figure 2 plus an exponential disk galaxy (see discussion in § 4). 

Radiation from the quasar will have a strong effect on the 
interstellar medium of the host galaxy (see, e.g., Begelman 
1985; Chang, Schiano, & Wolfe 1987). The radiation may 
strongly inhibit star formation near the quasar. Therefore, it is 
not obvious whether an aperture magnitude calculated 
assuming a constant surface brightness or an exponential disk 
model with a surface brightness increasing toward the quasar 
will give a more accurate estimate for the total magnitude of a 
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TABLE 5 


Estimated Characteristics of Candidate Host Galaxies 


Quasar 

m 606W 

{Aperture Magnitude) 

m 606W 

(Exponential disk) 

My 

Comments 

PG 1116 + 215 

18.6* 

18.1 

-20.1 (-20.7) 

Partial rings? 

3C 273 

17.1 b 

16.8 

-21.4 (-21.7) 

Off-center elliptical? 

PG 1444 + 407 

19. 4 C 

18.8 

— 20.3 ( — 20.8) 

Ring-like, bar? 


* Between 2" and 7". 
b Between 3" and 15". 
c Between 1"5 and 4". 


host galaxy. Moreover, our experiments with fitting model 
host galaxies to our quasar images suggest that the results are 
of uncertain significance. 

In the following discussion, whenever a best estimate for the 
magnitude of a host galaxy is required, we use the exponential 
disk magnitudes listed in Table 5. 

We also note (see Paper I) that for PG 0953 + 414 we see 
what may be an extended, low surface brightness feature (p v ~ 
24 mag arcsec 2 , between P.A. « 100° to ~P.A. % 180°) at 
~~ 3" to ~ 5" from the quasar image. 

6. COMPANION GALAXIES 

The unprocessed images of the quasar fields displayed in 
Figure 1 reveal a number of companion galaxies projected 
close to the quasars. Very close companions (separations <2") 
are shown more clearly in Figure 4, which presents the residual 
images after subtraction of a stellar PSF. 

In order to evaluate the statistical significance of the small 
projected separations of the quasar-galaxy pairs, we decided to 
count the number of companion galaxies brighter than a speci- 
fied limiting absolute magnitude that were found to have a 
metric separation from one of the quasars of less than or equal 
to some fixed distance. We choose a priori a limiting absolute 
magnitude of M v = —16.5 (4 mag fainter than L*) and a 
maximum separation of 25 kpc. All galaxies found around the 
quasars that satisfy these specifications are included in our 
complete sample of galaxy companions. 

Table 6 lists for each quasar the number of companion gal- 
axies that are at least as bright as M v = —16.5 and that are 
projected within 25 kpc of the center of light of the quasar. The 
table gives the separations both in arcseconds and, in paren- 
theses, in kiloparsecs. Similarly, the brightnesses of the galaxies 
are tabulated in both apparent and absolute magnitude. 

TABLE 6 

Galaxy Companions Brighter than M F606W = - 16.5 within 
25 kpc of the Quasar 


Distances Magnitudes 

Number of 

Quasar Companions arcsec kpc m M 


PG 0953+414 I 7.9 18.8 22.9 -16.5 

PG 1116 + 215 1 12.2 23.6 19.5 — 19.2 

PG 1202 + 281 2 5.0 9.2 19.3 -19.3 

9.3 17.1 21.9 -16.7 

3C 273 0 

PKS 1302-102 2 1.1 3.0 20.6 -19.2 

2.2 5.9 21.9 -17.8 

PG 1307+085 0 

PG 1444 + 407 0 ... 

3C 323.1 1 2.6 6.6 20.8 - 18.8 


How can we evaluate the statistical significance of the close 
pairs listed in Table 6? The ideal procedure would be to count 
galaxies in random areas that were also imaged with WFC2 
and the same filter used here, F606W. We intend to do this 
when the full sample of 20 quasar fields is available. In the 
meantime, we can establish an upper limit to the random prob- 
ability by counting galaxies in comparison regions imaged at 
the same time and with the same instrumental configuration as 
each of the quasar fields. We count galaxies in circular areas 
that are within 25 kpc (at the quasar redshift) of the quasar 
position shifted to chip 2 and chip 4. (As described earlier, see 
§ 2, the quasar images are all close to the center of chip 3 of the 
WFC2.) This procedure provides an upper limit to the random 
probability of finding galaxies of the specified minimum bright- 
ness within the given area because the centers of chip 2 and 
chip 4 are separated by only 80", about 170 kpc at z = 0.2, 
from the center of chip 3. The separation of the centers of chip 
2 and chip 3 (or chip 4 and chip 3) is, for objects at z = 0.2, 
rather similar to the typical core radius (~ 125 icpc; see Bahcall 
1977) of a rich Abell cluster of galaxies. Therefore, the compari- 
son regions in chip 2 and chip 4 may contain an enhanced 
number of galaxies if the quasar-galaxy correlation function is 
large at low redshifts. The areas inspected in the comparison 
regions are slightly larger than the accessible areas around the 
quasars, since companion galaxies are difficult to observe on 
our images if they are projected less than 1" from the quasar 
centers. 

In the 16 fields of chip 2 and chip 4 that accompany the eight 
quasar images studied here, there are six galaxy companions 
that would have satisfied the a priori criteria for inclusion in 
our complete sample of close companions had been found at 
the same coordinate positions in chip 3. 

We have also counted galaxies close to four other quasars 
that are in our complete sample of 20 quasars and for which 
data became available after the initial drafts of the present 
paper were finished. These new objects are PKS 2135—14 
(z = 0.200), PKS 2349-0.14 (z = 0.173), PHL 909 (z = 0.171), 
and NAB 0205 + 02 (z = 0.155). There are five additional close 
galaxy companions (two for PKS 2349-014 and one each for 
the other three new quasars) in the four new quasar fields; 
there are no new companion galaxies in the eight new compari- 
son regions (of chip 2 and chip 4). 

If we limit ourselves to the eight quasar fields studied in the 
present paper, then the average number of companion galaxies 
in chip 2 or chip 4 within the specified parameters is 6/16 or 
0.375. The Poisson probability of finding seven or more pro- 
jected close galaxy companions is listed in Table 6 is 

P(> 7; 8; 25 kpc, M < -16.5) < 0.02 . (6a) 

If we consider all twelve quasar fields imaged so far by HST , 
then the average number of companion galaxies in chip 2 or 
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chip 4 is 6/24 or 0.250. The Poisson probability of finding at 
least 12 projected companions observed close to the 12 quasars 
is 

P(> 12; 12; 25 kpc, M < - 16.5) < 0.0001 . (6b) 

We conclude that the presence of the projected close compan- 
ions seen in Figure 1 and Figure 4 is statistically significant. 

The very close companions, separations < 3”, are difficult to 
observe from the ground. There are three such very close com- 
panions of the eight quasars considered here that are brighter 
than our absolute magnitude limit of M y = —16.5 and none 
within 3" of the 16 equivalent quasar positions on chips 2 and 
3. If the three very close companions are not physically associ- 
ated with their neighboring quasars but instead are a chance 
superposition, then the post facto Poisson probability for 
finding no companions within 3" for the 16 comparison regions 
is 0.0025. There is one additional very close companion bright- 
er than M v = — 16.5 among the four new quasars (projected 
close to PKS 2349-014), but no galaxies within 3" of the 
equivalent quasar positions on chip 2 and chip 4. The post 
facto Poisson probability of finding a total of four very close 
companions (closer than 3") among the 12 quasar fields and 
zero very close companions in the 24 comparison regions is less 
than 0.01%, but post facto probabilities are often very small. 

The situation for PKS 1302-102 is especially remarkable 
and is shown in Figure 8 (Plate 20). Two companions of this 
quasar are seen clearly in Figure 8; they are projected at 
separations of only l'.T (3 kpc) and 272 (6 kpc) from the quasar 
location. They lie at projected distances that are well within 
what would be the volume enclosed by a typical L* host galaxy. 
Hutchings & Neff (1990) report the presence of an extremely 
luminous host galaxy, M y = —23.7, which we do not detect. 
Hutchings et al. (1994b) did not detect the companion at 272; it 
appeared to be an “ extended knot ” with the combined data 
from the original WFC and CFHT HRCam. 

7. BRIGHTNESS LIMITS ON MCLEOD-RIEKE HOST GALAXIES 

Recently, McLeod & Rieke (1994a, b) have performed a very 
careful ground-based study in the H band (1.6 fim) searching 
for host galaxies of quasars; all eight of the quasars discussed 
in this paper are in their sample. McLeod & Rieke (1994b) 
report detecting bright host galaxies in at least 23 of 26 high- 
luminosity, nearby quasars (M^< —23.1 mag and z < 0.3). 
They conclude that the high-luminosity quasars have relatively 
bright host galaxies, typically twice the luminosity of an L* 
(Schechter) galaxy. Their results are similar to those of a 
number of other workers who have used ground-based data in 
different wavelengths (see, e.g., Hutchings & Neff 1992; Veron- 
Cetty & Woltjer 1990 and other references cited in these papers 
and in McLeod & Rieke 1994a, b). 

Since the McLeod & Rieke observations were made in the H 
band and our observations are carried out in essentially the V 
band, we cannot compare our results directly with theirs. 
However, we can compute the average absolute magnitude in 
V to which we would be sensitive for host galaxies like those 
reported by McLeod & Rieke. By comparing our limit on 
(M y > with their reported <M H >, we can infer a minimum 
value of V — H for the two sets of results to be consistent. 

We have created simulated host galaxies which have the 
same exponential scale lengths as the host galaxies described 
by McLeod & Rieke (1994b) in order to determine our ability 
to detect hosts of the type they described. The simulated gal- 


axies we used are more difficult to detect than real galaxies 
with the same scale lengths. The simulated galaxies have no 
particular features that the computer can notice; they are 
smooth by construction. 

All of the host galaxies reported by McLeod & Rieke (1994b) 
for the quasars shown in Table 7 have apparent exponential 
scale lengths that are of order 1725 ± 0'.'5, i.e., typically of order 
2.5 kpc at the redshifts of the quasars. The seeing in their 
observations was also about 1'.'25 + 0.25 (see McLeod & 
Rieke 1994b). The similarity between the quoted seeing and the 
model scale lengths may be an indication that some caution 
might be required in estimating the systematic uncertainties. 
McLeod & Rieke did not publish their scale lengths because 
they used the exponential galaxy profiles as fitting functions, 
without claiming that the profiles represented real galaxies. We 
are grateful to K. McLeod for generously making available her 
unpublished values for the scale lengths. 

We created two simulated, face-on, smooth galaxies with 
exponential scale lengths of l'.'O and F5. We then performed 
tests identical to those described in § 3.3 to determine how 
faint such galaxies could be and escape detection in our tests. 

The results are shown in Table 7. On average, we are sensi- 
tive to the simulated host galaxies to V = 18.7, which corre- 
sponds, for the redshifts of the eight quasars we investigate, to 
M v = —20.4. Thus we are able to detect the smooth featureless 
simulated galaxies to 0.1 mag fainter than L*. Requiring our 
results to be consistent with those of McLeod & Rieke (1994b) 
would imply that the host galaxies of the eight quasars we 
studied are at least 1 mag redder in V — H than ordinary field 
galaxies. 

8. WHAT HAVE WE DONE WRONG? 

We are painfully aware that our initial results appear to be 
in conflict with suggested detections of relatively bright 
quasars by previous authors. Therefore, we have performed a 
number of tests of our analysis procedures to see if we have 
made any errors that might explain the difference between our 
results and previously published results. Unfortunately, we 
have not discovered an explanation. We list below some of the 
tests we did perform. 

We first compare in § 8.1 the ground-based PSFs used in 
previous work on host galaxies of quasars with the point- 
spread function measured for our HST observations. Next, we 
investigate in § 8.2 the possible uncertainties that result from 
errors in the sky subtraction and then determine in § 8.3 the 
likely uncertainties that are introduced by using a stellar point- 
spread function that might be somewhat different (because of 
spatial, temporal, or color changes in the PSF) than the PSF of 

TABLE 7 


Detectability of Smooth. Face-On 
Exponential Disk Galaxies 



Scale Length 

Quasar 

r:o 

175 

PG 0953 + 414 

18.5 

18.8 

PG 1116 + 215 

18.0 

18.7 

PG 1202 + 281 

19.1 

19.5 

3C 273 

17.0 

17.7 

PKS 1302-102 

18.5 

18.9 

PG 1307 + 085 

18.0 

18.7 

PG 1444-407 

18.9 

19.5 

3C 323.1 

19.2 

19.5 
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the nuclear region of the quasar. The intensity of the best- 
fitting stellar PSF may be determined in at least two extreme 
ways, either by minimizing the residual light by optimally 
matching the diffraction spikes or by azimuthally averaging 
the light in the quasar images. In § 8.4, we show that both 
methods yield essentially the same results. Next, we construct 
in § 8.5 pseudo-AGNs and show that we can detect the pres- 
ence of host galaxies approximately 4.5 mag fainter than the 
simulated quasars. Finally, in § 8.6 we examine compressed 
images and show that we have not missed large, low surface 
brightness features that might show up when the angular 
resolution was reduced. We also present general arguments 
showing that WFC2 images are capable of detecting low 
surface brightness objects that are much fainter than the (dark) 
HST sky level. 

8. 1 . Comparison of Ground-based and HST PSFs 

Can HST detect faint diffuse emission near a bright stellar 
source more effectively than previous ground-based observa- 
tions? The measurement of extended scattered light in the 
WFPC2 instrument (Krist & Burrows 1994) requires us to give 
a quantitative answer to this question. 

Figure 9 shows a direct comparison of the PSF produced by 
the WFPC2 (see Fig. 3) with the seeing profile in an excellent 
ground-based observation. Both curves present azimuthal 
averages. The HST PSF is shown in two parts; points inside of 
1" are from a theoretical model by Burrows (1994b), and points 
beyond 1" are from our composite PSF (Fig. 3). The ground- 
based seeing profile is from observations taken with the Hale 
telescope; the FWHM of the profile is 0787, and the half- 
energy diameter is 1701. 

At a radius of 1", the HST PSF has an advantage of about a 
factor of 10 over the ground-based Hale PSF. This advantage 
decreases somewhat at large radii because of scattering in the 
WFPC2 (Krist & Burrows 1994). 

How good is the seeing used in the previously published 
studies of the host galaxies of quasars? In the earliest papers on 
this subject (see references in § 1), the seeing conditions were 
often not reported, but in recent years most authors have 



Fig. 9. — Comparison of HST and ground-based (Palomar) point spread 
functions. This figure compares the azimuthally averaged radial surface bright- 
ness profile of the PSF of the F606W stellar HST images with the PSF 
obtained at the Hale Telescope with a V filter in 0"87 (FWHM) seeing. Both 
the (Burrows) phenomenological and the measured stellar PSF are shown for 
the HST data. 


TABLE 8 

Seeing Conditions (FWHM) Reported by Different 
Observers of Quasar Host Galaxies 


Observers 

Reported Seeing 
(FWHM) 

Wyckoff et al. 1981 

171-176 

Hutchings et al. 1982 

-1" 

Tyson et al. 1982 

175 

Malkan et al. 1984 

175-270 

Smith et al. 1986 

079-178 

Hutchings 1987 

079-173 

Hutchings et al. 1989 

077 -173 

Romanishen & Hintzen 1989 

171-274 

Veron-Cetty & Woltjer 1990 

079-172, 78% 
173-174, 22% 

Hutchings & Neff 1992 

-075 

Dunlop et al. 1993 

“ complex ” (076 pixels) 

McLeod & Rieke 1994b 

1725 ± 0725 


reported the FWHM of the seeing either quantitatively or 
qualitatively. Table 8 gives the FWHM of the PSF as reported 
in the previously published papers in which we found a definite 
statement about the seeing. The 0787 Hale Telescope shown in 
Figure 9 is better than nearly all the PSFs used in previous 
ground-based work. 

We conclude that the WFC2 PSF of the Hubble Space Tele- 
scope provides an important advantage over previous ground- 
based studies for projects in which one tries to detect faint 
diffuse emission near a bright stellar source. 

8.2. Sky Subtraction 

Oversubtracting the sky brightness would result in removing 
signal from a possible host galaxy. Sky oversubstraction is the 
astronomical equivalent of throwing out the baby with the 
bathwater. 

We have therefore measured the sky background over the 
fields of each of our CCD images. There are detectable (non- 
Poisson) variations in the flat-fielded sky intensity at the 
several percent level; the local deviations are smaller in the 
centers of the fields where quasar images are located than at 
the edges. (The sky levels in the corners of some fields are as 
much as 20% lower than at the centers of the fields.) 

The measured sky brightness for the eight quasar images we 
have studied range (probably due to different amounts of 
earthshine) from 22.2 mag arcsec" 2 (for 3C 273) to 23.0 mag 
arcsec 2 (for PG 1444-1-407). Given the measured small fluc- 
tuations in the sky brightness across a CCD in a given expo- 
sure, it is extremely unlikely that the estimated sky levels are 
incorrect by as much as 10%. Even a 10% error in the sky 
brightness would not affect significantly our inferred limits 
(Table 3) on the typical host galaxy magnitudes unless the host 
galaxies were as faint as m F606W = 21, which is much fainter 
than the reported magnitudes of host galaxies determined from 
ground-based measurements (see Table 7) and also fainter than 
our limiting sensitivity (see Table 3). 

We conclude that uncertainties in determining the sky 
brightness are probably not important in detecting host gal- 
axies at the magnitudes discussed in this paper. 

8.3. Subtracting One Stellar PSF from Another Stellar PSF 

How accurate is the PSF subtraction that we perform to 
detect the host galaxies? Could imperfect PSF subtraction 
cause us not to detect luminous host galaxies of the quasars? 
Are the residual diffuse features seen in Figure 4 associated 
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with the quasar or are they the result of imperfect stellar PSF 
subtraction? 

In order to answer these questions, we selected a representa- 
tive set of nine saturated stellar images that are contained in 
different WFC2 HST fields. The combined PSF of the stan- 
dard star F141 in M67 (see Fig. 2) was subtracted from the 
image of each one of the other nine stars. The subtracted 
images are displayed in Figure 10. 

The first six panels, a-e, of Figure 10 (Plate 21) show the 
results obtained when the standard PSF of Figure 2 was sub- 
tracted from one of the other saturated stellar images, both the 
standard PSF and the other saturated stellar image having 
been obtained with the F606W filter. The last three panels of 
Figure 10 show the residual images found when the standard 
PSF was subtracted from stellar images observed through the 
F814W, F658N, and F439W filters, respectively. 

The brightest residual diffuse emission is found in Figure 
10 d. This is also the case for which the original star had the 
largest saturated area of any of the test stars, considerably 
more than even in our longest exposures for the quasars. 
Moreover, the residual diffuse emission does not have the 
appearance of a galaxy. Among other things, the image con- 
tains a second, fainter star within 2" of the center of the bright 
star. The relatively poor subtraction shown in Figure 10 g also 
seems to be caused because the field star is highly saturated, 
again considerably more than our quasar images. 

The residual surface emission shown in Figure 10 is prob- 
ably an upper limit to what might be produced by inaccuracies 
in the stellar PSF subtraction from the actual quasar images. 
The reason it is an upper limit is that the standard stellar PSF, 
Figure 2, was obtained with the star observed through the 
same filter as the quasar and located at the same position as 
the quasar image; we purposely mismatched PSFs in order to 
produce Figure 10. 

We can estimate quantitatively the equivalent magnitudes of 
host galaxies that are produced by the mismatch between 
stellar PSFs, scaling the stars used to produce Figure 10 to the 
brightness levels that they would have if they were as bright as 
our average quasars, i.e., m F606 w = 15.2 mag. First, we calcu- 
late for each star the equivalent exposure time, f equiv , that a 
15.2 mag star would have been exposed in order to produce the 
number of photons counted for that star in the original image 
(before subtracting the standard PSF). Then we fit a model of 
the form given in equation (4) to the stellar light and determine 
the total number, JV host , of photons that are in the component 
assigned to the host galaxy. Finally, we calculate the magni- 
tude of the host galaxy using t equiv and N hosl . When exponential 
disk models were fitted to the HST stars, four of the pseudoga- 
laxies had positive counts and five had negative counts. The 
scale lengths were typically 20-40 pixels, and the pseudoga- 
laxies with positive flux had an average equivalent F606W 
magnitude of 

< m F606w( m i s match)> exp di9k = 21.2 , (7) 

and a median value of 20.8 mag. If one changes the sign of the 
negative flux pseudogalaxies and assigns a magnitude to this 
value, the average brightnesses of these five objects is 20.9. If 
we estimate instead all the light not in a point source using 
aperture magnitudes, we find considerably brighter residuals, 

<™ F60 6w(mismatch)> aper .mag 19.8 , (8) 


Even though the match between the standard and the field 
star PSFs is, by choice of the stellar images, not perfect, the 
subtraction is sufficiently accurate that there does not appear 
to be a significant danger that we have failed to observe lumi- 
nous host galaxies (m F606 w < 20.0) because of errors in the 
subtraction. 

8.4. Diffractions Spikes versus Azimuthal Average 

With ground-based observations, it is customary to deter- 
mine the intensity of the stellar image to be subtracted from the 
quasar image by azimuthally averaging the light intensity near 
to the quasar and then requiring that azimuthal average to 
match as well as possible the PSF for a point source. This 
procedure has the disadvantage that close companion galaxies 
(some of which are visible on HST images; see Figs. 1, 4, or 8) 
are not properly taken into account. We are therefore led to 
investigate whether the azimuthal average that is usually 
employed in ground-based investigations might be a major 
source of the discrepancy between our results and some of the 
previous studies. In this subsection, we answer the following 
limited questions: For the HST data, do we obtain signifi- 
cantly different estimates for the possible host galaxy magni- 
tudes if we use an azimuthally averaged radial surface 
brightness profile rather than concentrating on matching as 
well as possible the diffraction spikes? Our previous discussion 
in §§ 3 and 5 was based upon an approximately optimal 
matching of diffraction spikes in the point-source PSF and the 
quasar image. 

In order to answer the question posed here, we calculated in 
two ways the apparent magnitudes of host galaxies with expo- 
nential disks that, together with the composite point source 
PSF, best fit (in a x 2 sense) the quasar images. We determined 
the central intensity of the best-fit stellar PSF plus exponential 
disk by : (1) minimizing the residuals from the subtraction using 
all the emission in annuli centered on the quasar images; and 
(2) first minimizing the residuals from the quasar minus point 
source subtraction by considering only the diffraction spikes. 
The two methods gave consistent results: in no case was a 
candidate host galaxy seen when the subtraction was done one 
way but not seen when the subtraction was done the other 
way. In addition, the magnitudes found for the exponential 
disk galaxies were about the same using the two different tech- 
niques. On average, the magnitudes obtained when the diffrac- 
tion spikes were matched before the exponential disks were 
fitted to the residuals led to only slightly brighter “ host gal- 
axies.” Specifically, 

<m(azimuthal average) — m(spikes)) = (0.17 ±0.19) mag . (9) 

We conclude that the difference in subtraction procedures of 
concentrating on either the diffraction spikes or an azimuthal 
average does not cause us to overlook luminous host galaxies. 

8.5. Host Galaxies of Pseudo- AG Ns 

We have performed an overall test of our techniques for 
finding host galaxies by creating pseudo-AGNs within 
observed or simulated host galaxies and then applying our 
standard software to demonstrate that we can detect the host 
galaxies. The pseudo-AGNs were produced by adding host 
galaxies to a saturated image of a 15.5 mag bright star located 
13"8 from the center of chip 2 of the WFC2; this star was also 
used in producing Figure 10 The stellar PSF that was used for 
the AGN was chosen to be very different from the composite 
PSF shown in Figure 2 and, in fact, the image of the pseudo- 


and a median value of 20.4 mag. 
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AGN even has a small doughnut-shaped ghost image that is 
clearly visible in Figure 10 f We are therefore making it pur- 
posely difficult to detect the faint galaxy underneath the AGN 
by creating a poor match between the nuclear PSF and the 
composite standard PSF. The galaxy images that were used to 
make the simulated AGNs included HST-observed galaxies 
from Figure 5, galaxies observed with ground-based telescopes 
and adapted to HST conditions, and simulated (smooth, 
noiseless) galaxies. 

Each host galaxy was given a sequence of successively fainter 
magnitudes, usually beginning with 1. 5-2.0 mag fainter than 
the pseudonucleus. Just as described in § 3 for the analysis of 
the quasar images, we subtracted a best-fit standard PSF (Fig. 
2) for each case and then analyzed the results with the same 
software as was used in studying the quasars. We also inspect- 
ed visually the subtracted images, just as we did for the 
quasars. The results were similar to what was obtained in § 3 
for the quasar images plus the f/ST-observed galaxies from 
Figure 5. Moreover, the yf distribution as a function of the 
faintness of the host galaxy had the same form as for the real 
data (see Fig. 6). 

Figure 1 1 (Plate 22) shows the images obtained when a simu- 
lated, nearly edge-on exponential disk galaxy (exponential 
scale length of T.'5, approximately 3 kpc at 2 = 0.2) was added 
to the stellar source. The galaxy looks much like the HST 
galaxies we have observed except that it is smoother. In Figure 
11a, the galaxy is 1.7 mag fainter than the quasar nucleus (a 
15.5 mag stellar image). The galaxy is 0.5 mag fainter in each 
successive panel. The limit of detectability is 19.7 mag, which is 
the apparent magnitude of the galaxy in Figure 11 d. If the 
quasar nucleus in Figure 1 1 had the same absolute magnitude, 
M v — —23.9, as the average quasar in our sample, then we 
could detect the host galaxy of the simulated AGN down to 
approximately M v = — 19.7. 

We conclude that our technique is capable of discovering 
moderately bright to luminous host galaxies of quasars at 
small redshifts. 

8.6. Low Surface Brightness Galaxies 

Could we have missed low surface brightness galaxies that 
are luminous but which are spread out over a very large area? 
We have addressed this point earlier by analyzing simulations; 
we showed in § 3.3 that we could detect large simulated gal- 
axies to magnitudes comparable to what we were able to 
achieve for the more compact galaxies (see the last two rows of 
Table 3). In this subsection, we describe a direct test for the 
presence of low surface brightness images that were made 
using compressed images. We also present some general con- 
siderations which show that, with the moderately long expo- 
sures used in this paper, the WFC2 of HST is sensitive to 
relatively faint surface brightness objects. 

We formed compressed images from the original images, 
with or without the subtraction of a stellar PSF, by combining 
the light in 3 x 3 pixels into a single pixel. We then reexamined 
by eye all the images for the eight quasars analyzed in this 
paper. No features were observable in the compressed images 
that were not already visible in the uncompressed images. 
Except for 3C 273, the compressed images did not show a 
qualitative improvement in the visibility of faint diffuse fea- 
tures. The improvement in “seeing” the host galaxy in the 
compressed images of 3C 273 can be appreciated by comparing 
Figure 4 and Figure 7. For the quasar PKS 2349-014 (which 
we studied after the data discussed in the present paper were 


analyzed; see Bahcall et al. 1995), the large diffuse emission 
centered approximately 2'.'5 from the quasar center of light is 
also more clearly visible on the compressed images. However, 
for both these cases, the diffuse emission is unmistakably 
present in the uncompressed images. For PKS 1302—102, we 
are able to measure more easily on the compressed image a 
very faint plume (which lies almost north-south and is visible 
between about 15"-23" from the quasar) that has a very low 
surface brightness, about 25.2 mag arcsec' 2 , and a total appar- 
ent visual magnitude of about 19.1 mag (Af v = —20.7). 

The WFC2 images of HST can detect surface brightness 
levels well below the level of the sky. For an 1 100 s exposure 
taken with the filter F606W filter (the highest throughput 
filter), the sky level per pixel is about 125 electrons. The photon 
noise introduced by the sky is somewhat larger than the noise 
produced by the amplifiers in the CCD (seven electrons). This 
implies that one should be able to measure surface brightnesses 
down to a theoretical limit (set by the quantization of the 
analog-to-digital converters) of approximately g. v ~ 26 mag 
arcsec' 2 , although it is difficult to achieve this limit in practice 
because of uncertainties in the flat-field calibration and in 
other instrumental effects. We are able to see easily a large 
diffuse feature in the image of PKS 2349-014 that is [i v ~~ 24,5 
mag arcsec' 2 , about 2 mag below the average sky level. This 
situation is to be contrasted with observations taken with the 
planetary camera, WFPC2, whose image scale is such that it is 
always read-noise, not sky-noise, limited (each PC pixel con- 
tains only £ of the area of the WFC pixels). It is therefore very 
difficult to detect objects with low surface brightness with the 
PC. For detection of low surface brightness emission from 
quasar host galaxies, the PC2 is at a further disadvantage com- 
pared to the WFC2; the scattered light from a stellar source is, 
at the same angular distance, much higher in the PC2 than in 
the WFC2 (see Krist & Burrows 1994). 

9. SUMMARY AND DISCUSSION 

We have studied the environments of eight of the intrinsi- 
cally most luminous quasars {M y < -22.9) that are found at 
small redshifts (z < 0.3). The principal results of our analysis 
are summarized below. 

1. Separate companion galaxies are found projected so close 
to the quasar location that they fall well within the volume that 
would be occupied by an L* host galaxy (see § 6). For example, 
PKS 1302- 102 has two companion galaxies clearly visible in 
Figure 8 that lie at projected distances of 3 and 6 kpc, respec- 
tively, from the center of light of the quasar emission. Dynami- 
cal friction will cause these very close companions to spiral 
into the quasar in a time short compared to the Hubble time if 
the companions have typical galaxy masses. 

We have counted galaxy companions within a projected dis- 
tance of 25 kpc of the quasar position that are brighter than 
M v = — 16.5 for all eight of the quasars discussed in this paper 
and for four additional quasars for which data were obtained 
after this paper was almost complete. There are (see § 6) seven 
(12) companion galaxies found near the eight (12) quasars. By 
counting galaxies in comparison fields, we have established (see 
eqs. [6a] and [6b]) that most of the observed close galaxy 
companions are physically associated with the quasars. 

These statistical results provide strong evidence for the view 
(see references in § 1) that one important ingredient in the 
luminous quasar phenomenon is the presence of galaxy com- 
panions. The very close companions may be in the process of 
providing fresh fuel by spiraling into the quasar. 
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It will be of great interest, but difficult, to obtain redshifts for 
the companion galaxies that are listed in Table 6. The cosmo- 
logical nature of the quasar redshifts, assumed in this paper, 
plus the statistical association between quasars and compan- 
ion galaxies that is described above together imply that the 
many companion galaxy redshifts will be the same as the 
quasar redshifts. 

2. We have detected candidate host galaxies for three of the 
eight quasars we have studied, namely, for PG 1116 + 215, 3C 
273, and PG 1444 + 407. However, these candidate host gal- 
axies are not remarkably luminous; they have an average 
K-band brightness of about 0.6 mag brighter than L*. HST 
observations at a different telescope roll angles of these three 
quasars would provide crucial tests for the reality of the 
inferred characteristics, which are described in § 5, of the real 
candidate host galaxies. If the residual emission tentatively 
ascribed to host galaxies is real, it will have a fixed shape and 
position on the sky, independent of telescope orientation. 

Tyson et al. (1982) (see also Wyckoff et al. 1981; Roser & 
Meisenheimer 1991) used a 3" occulting disk to image 3C 273. 
They infer the presence of a luminous galaxy, V = 16.05 
{M Y = -22.4), with a major axis position angle of 70°, and a 
total extent of more than 40". The major axis of our candidate 
host galaxy for 3C 273 (M v = -21.7, see § 5 and Fig. 7) is 
consistent with the axis inferred by Tyson et al. Similarly, 
McLeod notes (McLeod 1995) that the faint halo reported by 
McLeod & Rieke (1994b) around PG 1444 + 407 has the same 
direction as the major axis (approximately north-south in Fig. 
1) of our HST candidate host galaxy. The agreement between 
the major-axis directions in the ground-based and the HST 
images constitutes important evidence for the reality of these 
candidate host galaxies. 

3. For the representative set of HST-o bserved galaxies 
shown in Figure 5, we could have detected host galaxies of the 
quasars if they had an average absolute magnitude as bright as 
M F606Vf (\im) = - 19.6 mag, i.e., to a limiting magnitude that is 
1 mag fainter than the characteristic (Schechter) magnitude L* 
for field galaxies. This average limiting magnitude was estab- 
lished in § 3.3 by placing a sequence of progressively fainter 
versions of each of the galaxies shown in Figure 5 underneath 
each of the quasar images and analyzing the residual emission 
that is left over after a point-source PSF was subtracted from 
the quasar plus added galaxy. 

4. We have fitted smooth, featureless galaxy models 
(exponential disks or de Vaucouleurs profiles) to the residual 
light after a best-fitting stellar PSF was removed from the 
quasar images. The fits are generally poor, but they do reveal 
(§ 4) evidence for possible host galaxy light that, on the 
average, is about equal to L*. This is also about the limit of 
detectability found in § 3.3 for early-type HST-observed gal- 
axies. 

5. Three of the luminous nearby quasars we have studied 
have apparent host galaxies that are moderately luminous and 
the other five quasars have no obvious bright host galaxies. In 
addition, PKS 2349-014, another member of our luminous 
quasar sample, has (see Bahcall et al. 1995) a very complex 
environment which includes apparent tidal wisps and a huge, 
off-center (from the quasar nucleus) low surface brightness 
nebulosity. Some of the quasars in our sample have apparent 
galactic companions projected very close to the quasar center 
of light (see Table 6). We conclude that luminous nearby 
quasars exist in a variety of environments. 

6. Our results on the absence of luminous host galaxies 


around some of the luminous quasars are surprising given the 
generally accepted view that quasars reside in host galaxies 
with the most luminous quasars residing in the most luminous 
galaxies. It is difficult to compare our results directly with 
previous ground-based observations because of the higher 
angular resolution available in the HST images and, in many 
cases, because of the different photometric bands (or spectra) 
that were used. It is not clear how the presence of close com- 
panions (revealed by HST images) would affect measurements 
from the gound of the azimuthally averaged residual light dis- 
tributions remaining after a point-source PSF is subtracted. 

To illustrate some of the considerations that should be taken 
into account in a comprehensive study, we have compared the 
HST results with several ground-based observations for which 
previous observers have published quantitative results for 
quasars in the sample studied in the present paper. For 
example, we have made a global comparison of our analysis 
with the recent results of McLeod & Rieke (1994b), who have 
reported finding bright host galaxies in H-band images for all 
eight of the quasars studied in the present paper. In order for 
our results to be consistent with the McLeod-Rieke conclu- 
sions, the host galaxies of luminous quasars must, on the 
average, be at least 1 mag redder in V - H than normal field 
galaxies. Hutchings et at. (1989), on the other hand, have 
reported that the host galaxy of PG 0953 + 414 has M B = 
-20.9, which would be consistent with our upper limit if the 
host galaxy has an extremely blue color, i.e., B-V < —1.0 
mag. 

Hutchings & Neff (1990) report, based upon ground-based 
observations, the presence of an exceptionally luminous host 
galaxy for PKS 1302- 102; they give an absolute magnitude of 
M = - 23.7, much brighter than the brightest galaxies of rich 
clusters. Although their observations were taken in R and ours 
were taken in V , we have not been able to think of a plausible 
explanation of why our upper limits for the host galaxy light 
are so much less than their detection. Our average upper limit 
for host galaxies like the HST-observed galaxies in Figure 5 is 
My = —20.1 (see Table 1), and the exponential disk fit to the 
total residual light is M v = — 21.4 (see Table 4). 

Veron-Cetty & Woltjer (1990) have suggested that the 
apparent magnitudes of host galaxies should be measured in a 
fixed metric annulus that is well removed from the quasar 
nucleus. They propose an annular region of 12.5-25.0 kpc for 
Q 0 = 0.0 and H 0 = 50 km s” 1 Mpc -1 . We have only one 
object in common with Veron-Cetty & Woltjer, PKS 
1302-102. The annular region specified by Veron-Cetty & 
Woltjer is between 2'.'175 and 4735 for PKS 1302— 102, which 
includes contributions from the two faint companion galaxies 
shown in Figure 8. The F-magnitude that Veron-Cetty & 
Woltjer (1990) obtain after transforming their /-magnitude 
measurement for this annulus is V — (19.46 + 0.2) mag. We 
obtain m F606 w = 19.6 ± °- 2 > in good agreement with their 
result. However, Veron-Cetty & Woltjer (1990) estimate that 
the host galaxy has a total M v = —23.7 for Q 0 — 0.0 and 
H 0 = 50 km s" 1 Mpc" 1 , which corresponds to M v = —22.0 
for the cosmological parameters used in this paper. The host 
galaxy reported by Veron-Cetty & Woltjer for PKS 1302— 102 
is as bright as the brightest galaxy of a rich cluster. For the 
same quasar, we are only able to set an average upper limit to 
the brightness of a host galaxy of M y = —20.1 (see Table 1) if 
the hosts are like the HST sample shown in Figure 5. On the 
other hand, we obtain host galaxy luminosities in the range 
M v = - 21 to - 22 if we fit the residual light to an exponential 
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disk or a de Vaucouleurs profile and try various initial guesses 
for the model parameters. Because of difficulties in interpreting 
the results of fits to smooth model galaxies (see § 4), and 
because of the presence of close companion galaxies revealed 
by HST (see Fig. 8), we are unable to reach any definitive 
conclusion in this case. 

7. Most of the work in preparing this paper was devoted to 
trying to answer the question : What have we done wrong? The 
longest section in this paper, § 8, describes a series of tests of 
our methods of analysis. These tests include: (1) measurements 
of the systematic errors due to sky subtraction ; (2) a determi- 
nation of the likely errors due to using a stellar PSF that is 
somewhat different from the PSF of the AGN nucleus; (3) 
adjusting the central intensity of the stellar PSF to either the 
diffraction spikes or the light within an annular region; and (4) 
constructing and analyzing simulated AGNs whose properties 
are known. Our analysis procedures successfully pass all these 
tests. Some of the techniques we have described in § 8 for 
testing the validity of our methods of analysis may also be 
useful in connection with ground-based observations. 

8. Although our results are not consistent with the most 
luminous quasars residing in the most luminous known gal- 
axies, our results are consistent with luminous quasars 
occurring in galaxies with a Schechter luminosity function that 
is cut off at a moderate luminosity (Fall 1994). We have evalu- 
ated the probability that we would have obtained the eight 
brightness upper limits given in the last column of Table 1 (or 


more stringent limits), if the host galaxies have a Schechter 
luminosity function (see Efstathiou et al. 1988) that is cut off on 
the low-luminosity side of 0.1L* (or 0.0 1L*). The probability of 
obtaining a result as stringent as we have found is somewhat 
unlikely, 0.06, if the lower cutoff in the luminosity function is 
0.1L* (average host galaxy luminosity of 0.5L*) but is as large as 
0.39 if the lower cutoff is 0.0 1L* (average host galaxy lumi- 
nosity of 0.2L*). The probabilities are even larger if we use the 
weaker limits that are obtained in § 4 for smooth featureless 
galaxies. In this case, we can approximate all the upper limits 
by L*. The probability of finding eight upper limits equal to L* 
(or fainter) is 0.4 (0.7) if the lower cutoff is 0. 1 (0.0 1 ). 
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Fig. 1. The unprocessed images of eight luminous quasars. This figure shows the long-exposure images of the eight luminous quasars observed with WFC2 
through filter F606W. 
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Fig. 4. The eight quasars with a best-tit stellar PSF subtracted. This figure shows the '‘long” (1400 s for PC i 1307 + 085; 1100 s for the others) F60bW 
observation of each of the eight quasars discussed in this paper. A best-fit stellar image from the Mb? calibration data has been subtracted from each quasar image. 
Each of the panels is 20" x 20"; the image scale is 070966 pixel 1 Note that 1 " ^ 2 kpc (0.2 r). 
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Fio. 4 — Continued 
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Fig. 5 . — HST images of eight galaxies. This figure shows the eight galaxies, taken from the fields of the quasar WFC2 observations, that were used to determine 
the sensitivity of the HST images to different types of quasar host galaxies. Each panel is 1 2'.'5 on a side. 
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Fig. 11 —Simulated AGNs. The figure shows a simulated AGN created from a bright stellar PSF (for the AGN nucleus) and a simulated nearly edge-on disk 
galaxy with an exponential scale height of 175. The stellar PSF is taken from a slightly saturated star located 14" from the center of CCD 2. The pseudo host galaxy 
has an apparent magnitude of 18.2 in the upper left panel and an apparent magnitude of 19.7 in the lower right panel. The limit of detectability with our techniques is 
19.7 mag. 
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